Moment estimates for solutions of linear stochastic 
differential equations driven by analytic fractional 

Brownian motion 



Jeremie Unterberger 



As a general rule, differential equations driven by a multi-dimensional irreg- 
ular path r are solved by constructing a rough path over F. The domain of 
definition - and also estimates - of the solutions depend on upper bounds 
for the rough path; these general, deterministic estimates are too crude to 
apply e.g. to the solutions of stochastic differential equations with linear 
coefficients driven by a Gaussian process with Holder regularity a < 1/2. 

We prove here (by showing convergence of Chen's series) that linear stochas- 
tic differential equations driven by analytic fractional Brownian motion [7l|8] 
with arbitrary Hurst index a £ (0, 1) may be solved on the closed upper half- 
plane, and that the solutions have finite variance. 
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Introduction 

Assume = (rt(l), . . . , Tt{d)) is a smooth d-dimensional path, and Vi, . . . , : 
M*" — > W be smooth vector fields. Then (by the classical Cauchy-Lipschitz 
theorem for instance) the differential equation driven by T 

d 

dy{t) = Y,V^{y{t))dTi{t) (0.1) 

i=l 

admits a unique solution with initial condition y(0) = yo- The usual way to 
prove this is by showing (by a functional fixed-point theorem) that iterated 
integrals 

{t) ■=yo+ } Viiyn{s))dTiis) (0.2) 
Jo . 

converge when n — > oo. 

Expanding this expression to all orders yields formally for an arbitrary 
analytic function / 

oo 

fiVt) = fiVs) + Y1 E [^^1 • • • ^^"/] iys)mh, in), (0.3) 

n=l l<ii,...,iri<d 

where 

/* rh rtn~i 

dTtAh) / dTt,{i2)... / dTtAin), (0.4) 

provided, of course, the series converges. By specializing to the identity 
function / = Id : M*" — > M'~, a; — > x, one gets a series expansion for the 
solution (yt). 

This formula, somewhat generalized, has been used in a variety of con- 
texts: 

1. Let 

N 

£^^'''{ys)=ys + Y. E [Vn■..V,Jd]{ys)m^,,...,^n) (0.5) 

n=l l<ii,...,in<d 

be the A^-th order truncation of (j0.3p . It may be interpreted as one 
iteration of the numerical Euler scheme of order A^, which is defined 

by 

2/r"'^'':=<'*"*'^-^o...of?'*^'*°(yo) (0.6) 



2 



for an arbitrary partition D = {0 = tQ < . . . < tn = T} oi the interval 
[0, T] . When F is only a-Holder with < a < ;^ , the iterated 

integrals T"'{ii, . . . , in), n = 2, . . . , N do not make sense a priori and 
must be substituted with a geometric rough path over T. A geometric 
rough path over F is a family 

{i^ts{h))i<h<d, (i^'tsik,i2))i<ii,i2<d, ■■■ , (Xtiih, ■ ■ ■ ,'i'N)i<ii,...,iN<d)) 

(0.7) 

of functions of two variables such that: 

(i) ri =Tl-Tl; 

(ii) {Holder continuity) each component of F^, A; = 1, . . . , is ka- 

Holder continuous, that is to say, sup^^j^ ^sup^gjj ^-^j|r^pr^^ < 
oo. 

(iii) [multiplicativity) letting ST^^^ := Fj^ — F^„ — F^^, one requires 

5FL(n,...,ifc)= E ^tik,...,ik,)riUik,+i,...,ik). 

ki+k2=k 

(0.8) 

(iii) {geometricity) 

T^^{il, . . . ,ini)T^^{ji, . . . jjn^) = F"i'^"2(/ci, . . . , fc^^+nj) 

fcesh{i,j) 

(0.9) 

where Sh{i,j) is the subset of permutations of ii , . . . , i„j , ji , . . . , 
which do not change the orderings of (ii , . . . , i„ J and (ji , . . . , jn^ ) ■ 

Properties (i)-(iv) are true when F is regular; the multiplicative prop- 
erty measures in some sense the defect of additivity of iterated inte- 
grals, which is easy to measure when one represents these as geometric 
quantities (areas, volumes, etc.) Under these conditions, it is possible 
to integrate a 1-form along the path F (or, more precisely, along the 
rough path F); we refer the reader either to [2j or to [3]. 

It is also possible to solve differential equations driven by F like (10. ip . 
either by using eq. (|0.2p and a fixed-point theorem in a class of 
F-controlled processes [3j, or by using the above Euler scheme [2j. 
Assuming either (i) the vector fields V and their derivatives up to 
order are bounded or (ii) they are linear 0, then the solution is 

^Similar results are also expected when the vector fields V are uniformly Lipschitz on 
R [4]. See also [5] for investigations on the same subject using a fixed point method. 
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globally defined, and the solution at time T is bounded (i) by a 
polynomial in |||r||| or (ii) by something like expC|||r|||^, where 

llirill = max„=i^..._Arsupi<j^^...^j^<rfSupo<^^t<T is the max- 

imum of the Holder norms. It seems there is no way (using either 
approach) to improve these bounds in the general deterministic set- 
ting. Unfortunately, they do not give a control of the solution as a 
stochastic process in the linear case (ii) when F is a Gaussian pro- 
cess (such as fBm or analytic fBm, see below) with Holder regularity 
a < ^, since EexpC(rt — F^)^ = 00 for C large enough, and in any 
case EexpClTf - T^l^ = 00 if iV > 3. 

2. Assume F is a stochastic process, and let Pt{f) = Kf{yt). When (jO.ip 
is a diffusion driven by usual Brownian motion, Pt is the associated 
semi-group operator. Assume now F is more general, for instance 
£Bm or analytic fBm. Even though the process in not Markov, the 
operator Pt is interesting in itself. The small-time expansion of Pt 
(corresponding to an arbitrary truncation of the above series) has been 
studied [1] when F is fBm with Hurst index a > 1/3. When a > 1/2, 
it has been proved [6] that the series converges for functions and vector 
fields V satisfying somewhat drastic conditions. 

In any case, it seems difficult to get moment estimates for the solutions 
of stochastic differential equations driven by stochastic processes F with 
Holder regularity a < ^. One of the reasons [S] is the difficulty of getting 
estimates for the iterated integrals F; another reason lies in the essence of 
the rough path method which relies on pathwise estimates; a third reason 
is, of course, that the Chen series diverges even in the simplest cases (one- 
dimensional usual Brownian motion for instance) as soon as the vector fields 
are unbounded and non-linear, e.g. quadratic. 

In this article, we prove convergence of the series (j0.3p when the vector 
fields Vi are linear and F is analytic fBm (afBm for short). This process - 
first defined in [8] -, depending on an index a £ (0, 1), is a complex- valued 
process, a.s. K-Holder for every k < a, which has an analytic continuation 
to the upper half-plane H"*" := {z = x + iy \ x £ 'R,y > 0}. Its real part 
(2Re Fj,i E M) has the same law as fBm with Hurst index a. Trajectories 
of F on the closed upper half-plane II"*' = H^ U M have the same regularity 
as those of fBm, namely, they are K-Holder for every k < a. As shown in 
[7], the regularized rough path - constructed by moving inside the upper 
half-plane through an imaginary translation i — > t + ie - converges in the 
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limit e ^ to a geometric rough path over T for any a £ (0,1), which 
makes it possible to produce strong, local pathwise solutions of stochastic 
differential equations driven by T with analytic coefficients. 

We do not enquire about the convergence of the series (j0.3p in the general 
case (as mentioned before, it diverges e.g. when V is quadratic), but only 
in the linear case. One obtains: 

Main theorem. 

Let Vi, . . . ,V(i be linear vector fields on C. Then the series \0.3\) con- 
verges in L'^{Q) on the closed upper half-plane n"*" = UM. Furthermore, 
the solution {yt)t<^u+' defined as the limit of the series, has finite variance. 
More precisely, there exists a constant C such that 

IE|yt-ys|' < C|t-s|2°, s,tGU+. (0.10) 

Notation. Constants (possibly depending on a) are generally denoted 
by C, C", Ci,Ca and so on. 

1 Definition of afBm and first estimates 

We briefly recall to begin with the definition of the analytic fractional Brow- 
nian motion F, which is a complex- valued process defined on the closed upper 
half-plane 11"'" [7]. Its introduction was initially motivated by the possibility 
to construct quite easily iterated integrals of F by a contour deformation. 
Alternatively, its Fourier transform is supported on M+, which makes the 
regularization procedure in [9l |10] void. 

Proposition 1.1 Let {^1,^1; A: > 0} be two families of independent stan- 
dard Gaussian random variables, defined on a complete probability space 
{Q,T,¥), and for k > 0, set = £.1 + i^fc- Consider the process F' defined 
for z eU+ byT'^ = Ylk>o fk{z)^k- Then: 

(1) F' is a well-defined analytic process on , with Hermitian covariance 
kernel 

EF'X = 0, EFX = 4^^1^(-i(z-«)))2"-2. (1.1) 

2 cos na 

(2) Let 7 : (0, 1) — > n"*" be any continuous path with endpoints 7(0) = 
and 7(1) = z, and set F^ = j^T'^du. Then F is an analytic process on 11"^. 
Furthermore, as z runs along any path in 11^ going to t £ M, the random 
variables F^ converge almost surely to a random variable called again Fj. 
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(3) The family {Fj; t G M} defines a Gaussian centered complex-valued pro- 
cess, whose covariance function is given by: 

_ p—'nra sgn(s) I c? 1 2o i pina sgn(t) I + 1 2a piTro sgn{t— s) I „ f 1 2a 

E[r,rt] = o, E[Tsft] = ^ — ^- ^ 

4cos(7raj 

The paths of this process are almost surely K-Holder for any k < a. 

(4) Both real and imaginary parts of {Tf] t G R} are (non independent) 
fractional Brownian motions indexed by R, with covariance given by 

E[Re rjm Ft] = J-^^ [-sgn(s)|s|2° + sgn(t)|t|2- - sgn(t - s)\t - s^"] . 

(1.2) 

Definition 1.2 LetYt := Re Tu, t G R+. More generally, Yt = (Yt{l), ■ ■ ■ ,Yt{d)) 
is a vector-valued process with d independent, identically distributed compo- 
nents. 

The above results imply that Yt is real-analytic on M^. 
Lemma 1.3 The infinitesimal covariance function ofYt is: 

EyX = 4^^(. + t)2"-2. (1.3) 

4 cos vra 

Proof. Let Xt := Im Tu. Since EF^Lt = 0, {Ys,s > 0) and {Xs,s > 0) 
have same law, with covariance kernel Wy^Yt = EX^Xt = iRe TisTu. Hence 

E[y;y/] = iRe EL^f = ^ii:M(, + tf^-\ (1.4) 
z 4 cos vra 

□ 

Note that WY^Yl > 0. From this simple remark follows (see proof of a 
similar statement in [6] concerning usual fractional Brownian motion with 
Hurst index a > 1/2): 

Lemma 1.4 Let YJ^(ii, . . . ,i„); n > 2 be the iterated integrals ofY. Then 
there exists a constant C > such that 

VarYr,(ii, ...,in)< (1.5) 

n! 
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Proof. Let 11 be the set of all pairings tt of the set {1, . . . , 2n} such that 
{{ki, ^2) G tt) =^ ^ifcj = ^fc^^ ; where k[ = ki if ki < n, ki — n otherwise, and 
similarly for fcg. By Wick's formula, 

VarYi",(ii,...,i„) 

/ dxi... dXn j i / dXn+l ... I dx2n 

^ S 'J S / \"/ S s 

n (1.6) 

Since the process Y' is positively correlated, and 11 is largest when all indices 
zi, . . . , are equal, one gets VarYj^(ii, . . . ,in ) < VarYi"^(l,...,l). On the 
other hand, Yf,(l, ...,l) = ^{Yt- Y,)'", hence 

V.rY«n ,._ [y^<Y,-Y,)r (2n)! ^ [2Var(y, - Y,)]" 
VarY,,(l, . . . , 1) - . < . (1.7) 

Now (assuming for instance < s < t) 

YaT{Yt -Ys)=Ca [ f {U + vf"-'^dudv < Cas'^''-^{t - sf < Ca{t - s) 
J s J s 



if I < s < t, and 



2a 

(1.8) 



2a[2a — 1) 

(1.9) 

if s < t/2. Hence the result. □ 

2 Estimates for iterated integrals of F 

The main tool for the study of T is the use of contour deformation. Iterated 
integrals of T are particular cases of analytic iterated integrals, see [8] or [7]. 
In particular, the following holds: 

Lemma 2.1 Let 7 : (0, 1) 11+ be the piecewise linear contour with affine 
parametrization defined by : 

(i) 7([0,l/3]) = [s,s + i|Re {t-s)\]; 

(11) 7([l/3,2/3]) = [s + i|Re (t - s)|, t + i|Re {t - s)\]; 
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(m) 7([2A1]) = [t + i\Re {t - s)\,t]. 

If z = 7(x) G 7([0,1]), we let 7^ be the same path stopped at z, i.e. 
Iz = 7([0, a;]), with the same parametrization. Then (letting Ca = 2ro'sTO ^ 

Varr[',(ii,...,i„) = 

c« X] / / ^"'iC-K^i - '^^(1)))^""^ • / dz2 dw2i-i{z2 - w^^2))) 



2a-2 



(iz„ / di(;„(-i(zn - 'u;a(„))) " (2.1) 

where S/ is i/ie subset of permutations of {!,..., n} such that {ij = 
ik) {(tU) =(7{k)). 



Proof. Note first that, similarly to eq. (jl.6 

rt rz-n-i \ f rt f-Wn-l 



/ dzi... dznji dwi ... I dWn 

ut^/ ^ ^ \ ^ ^ / 

n 

(2-2) 

(the difference with respect to eq. ()1.6p comes from the fact that contractions 
only operate between F's and F's, since E[r2^,r^j.] = Ep^u^T^u^,] = by 
Proposition II. ip . Now the result comes from a deformation of contour, see 

m- □ 

Lemma 2.2 There exists a constant C such that, for every s,t G n+ = 

n+ uM, 

Varrr,(n, ...,in)< ^ ' / • (2.3) 

n! 

Proof. We assume (without loss of generality) that Im s < Im t. If 
|Im (t — s)| > cRe \t — s\ for some positive constant c (or equivalently |Re (t — 
•s)| < c'\t — s\ for some < c' < 1) then it is preferable to integrate along 
the straight line [s,t] = {z £ C \ z = {1 — u)s + ut, < n < 1} instead of 
7, and use the parametrization y = Im z. U zi, Z2 G [s,t], yi = Im Zi,y2 = 
Im Z2, then |(-i(zi - ^2))^""^! < C(yl + y2)^""^ hence Varr^*,(n, i„) < 
C^VarY^^,^! (^1; • • • )^n)) which yields the result by Lemma [L^ So we shall 
assume that |Re (t — s)| > c\t — s\ for some constant c > 0. 
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Let us use as new variable the parametrization coordinate x along 7. 
Then formula (j2.ip reads 



VavTl^iii,. . . ,in) = c"^ y2 / '^^i / ^^i-^'C^i' 2/^7(1)) • / dx2 dy2K'{x2,ya{2), 

JO Jo Jo Jo 



dXn / dy-aK { 

-^n-i Vain )), (2.4) 

Jo 

where /s:'(x,2/) = (3|Re {t - s)|)2(3(x + y)|Re {t - s)\ + 2Im s)2"-2 q < 
x,y < 1/3, (3|Re - s)\f{3{{l - x) + {I - y))\Re {t - s)\ + 2Im t)2"-2 
if 2/3 < x,y < 1, and is bounded by a constant times \t — sp" otherwise 
thanks to the condition |Re {t — s)\ > c\t — s\. Note that (x + > 
22a-2 if < x,y < 1. Hence (if < < 1) \K'{x,y)\ < (Ci|t - 
[{x + + ((1 - x) + (1 - , which is (up to a coefficient) 

the infinitesimal covariance of \t - s|°(y^ + < X < 1) if y = ' y 

is independent of Y. A slight modification of the argument of Lemma 11.41 
yields 



Varrr,(n,...,i„) < {Ci\t-s 



^2„, [Var(yi - Yo) + Var(yi - YoT (2n)! 



(n!)2 ■ 2*^ . n! 

, (2CH-.|)^"°^ 

n! 

□ 



3 Proof of main theorem 

We now prove the theorem stated in the introduction, which is really a 
simple corollary of Lemma 12.21 

Let C be the maximum of the matrix norms |||Vi||| = sup||^||^^]^ ||Vix||oo 
for the supremum norm ||x||oo = sup(|xi|, . . . , \xr\)- Rewrite eq. (10. 50 as 

N 

£^^''\ys) =ys + Y^ «n,...,i.rr,(ii, . . . ,i„). (3.1) 

n=l l<ii,...,i„<(i 

Then ||ai,„„,ij|oo < and E|rj^,(ii, . . . , i„)|2 < (^^z^. Hence (by 
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the Cauchy-Schwarz inequality) 



AT 



{C"\t - S I 



E 




m,n=l 



( 



N 



) 



2 



.m=l 



E 




ma 



< C"".|t-s|2° (3.2) 



independently of A''. The series obviously converges and yields eq. (j0.10|) 
for p = 1. □ 
It should be easy to prove along the same lines that the series defining 
¥,\yt — y^p^ converges for every p > 1, and that there exists a constant Cp 
such that E\yt - y^pP < Cp\t - sp"^ for every s, t e n+. 
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